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Visual motion cues are one of the most important factors for eliciting animal behaviour, including 
predator-prey interactions in aquatic environments. To understand the elements of motion that cause such 
selective predation behaviour, we used a virtual plankton system where the predation behaviour in response 
to computer-generated prey was analysed. First, we performed motion analysis of zooplankton {Daphnia 
magna) to extract mathematical functions for biologically relevant motions of prey. Next, virtual prey 
models were programmed on a computer and presented to medaka (Oryzias latipes), which served as 
predatory fish. Medaka exhibited predation behaviour against several characteristic virtual plankton 
movements, particularly against a swimming pattern that could be characterised as pink noise motion. 
Analysing prey-predator interactions via pink noise motion will be an interesting research field in the future. 

For many predatory animals, survival depends on their ability to properly and rapidly identify objects, such as 
high-value prey. While hunting, a predator has to distinguish living prey from other environmental objects 
in real time. Naturally, all sensory information is used in this predator-prey interaction, including visual 
cues, which are an especially important factor in many cases 16 . In aquatic ecosystems, for example, selective visual 
predation by planktivorous fish greatly affects the composition of the zooplankton community 7,8 . Planktivorous 
fish select prey based upon numerous visual factors, including the size, visibility, colour, shape, and motion of the 
plankton 913 . Several previous studies have postulated that the visual motion of zooplankton is the most important 
for inducing predatory behaviour 814-16 . 

Only a few studies have examined predation behaviour based on prey swimming patterns; this is mainly due to 
the difficulty of controlling zooplankton swimming behaviour. In 1995 and 2005, a virtual plankton system, in 
which computer-generated prey was used to investigate vision in predator animals, overcame such problems and 
clarified the visual motion patterns affecting the induction of predation behaviour 14,15 . Bluegill sunfish {Lepomis 
macrochirus) preferred virtual plankton displaying higher speeds 14 or particular temporal motion patterns, such 
as a spinning motion 15 . However, mathematical modelling for the underlying spatiotemporal structure of the 
entire prey motion has not yet been attempted. Therefore, we have statistically analysed the swimming patterns of 
zooplankton (the water flea Daphnia magna) and determined the validity of the modelling. The predation 
behaviour of fish (the medaka Oryzias latipes) against computer-generated virtual plankton was studied. 

Medaka are known to have highly acute vision, as demonstrated by previous studies on the optomotor 
response 17,18 and open-field testing 1 '. Medaka are omnivorous fish, but dissection studies of their alimentary 
tracts have revealed that the gastric contents of medaka living in a natural pool contained up to 30% zooplankton 
and small aquatic insects, suggesting that medaka act as predators of self-moving small animals 20 . In the fresh- 
water environment of Japan, small Crustacea, including Daphnia, comprise the largest zooplankton population 
after Rotatoria 21 , and they play an important role in nourishing invertebrate and vertebrate predators, including 
medaka. Therefore we used medaka and water fleas {D. magna) as model organisms to study predator-prey 
interactions. 

The majority of zooplankton is known to swim in a random fashion 22 with irregular changes in speed and 
trajectory. This randomly complicated time-series of movements can be analysed by Fourier analysis, in which 
different wave sources are transformed to different characteristic statistical properties. These properties are called 
power spectrum densities (PSDs, power distributions in the frequency spectrum) and can be used to distinguish 
different types of waves. This classification by PSD is named using "colour" terminology, such as white noise 
(a signal with a flat frequency spectrum), pink noise (a density proportional to the inverse of its frequency) and 
blue noise (a density that is proportional to its frequency). Colour noises are widely found in nature. Pink noise, 
for example, has been found in neuron dynamics 23 . Retinal cells are arranged in a blue noise-like pattern, which 
yields good visual resolution 24 . In the present study, zooplankton swimming was classified using Fourier analysis. 
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Computer models of zooplanktons were then generated using an 
inverse Fourier transformation. Our data suggest an important role 
for the spatiotemporal motion patterns that are characterised as pink 
noise. 

Results 

Motion analysis of D. Magna. Figure 1A shows the swimming 
pattern of 1 plankton that travelled 334.40 mm from its start 
position to its end position over a 57 s period (1700 frames, mean 
speed = 5.87 mm/s). The swimming pattern of D. magna appeared 
to be a type of irregular walk; the moving distance and direction per 
frame changed constantly. Occasionally, D. magna suddenly stopped 
all motion, behaved like a dead particle, and appeared to sink due to 
gravity. 

To quantify the swimming pattern of the plankton, coordinate 
data from the first 257 frames (256 vectors of swimming velocity) 
of observation were taken for 48 plankton and summarised in Fig. IB 
as a frequency distribution of the swimming velocity (12,288 vectors 
in total). The results clearly illustrate that the frequency distribution 
has a single peak with wide tails located at either side. However, an 
asymmetric frequency distribution pattern was found in the vertical 
direction. Although the distribution frequencies gradually decreased 
as the speed increased in the upward direction, the distribution fre- 
quencies in the downward direction were concentrated at low speeds 
and formed a single peak around —8.0 mm/s. In contrast, a sym- 
metric frequency distribution pattern was found in the horizontal 
direction. The distribution frequency gradually decreases as the 
speed increases in both the left and right directions. These conclu- 
sions were also supported by the mean speed in each direction 
(Fig 1C, left graph). An approximately 2-fold difference in speed 
was found between the mean speeds in the downward and upward 
directions (downward, 4.41 ± 0.28 mm/s; upward, 7.67 ± 0.30; p < 
0.01, unpaired f-test). A smaller difference was observed between the 
mean speeds in the left and right directions (left, 3.95 ± 0.27 mm/s; 
right, 4.92 ± 0.29; p < 0.01, unpaired f-test). A similar result was also 
found for freshwater copepods {Cyclops abyssorum praealpinus) 15 . 
This asymmetric zooplankton motion pattern, including the one 
shown here for daphnids, reflects a bias to sink rather than to swim 
downwards 25 . A pair of swimmerets (the second antennae of D. 
magna) generates a preferential driving force in the upward dir- 
ection, and as a result, the main cause of downward motion must 
be a passive fall due to gravity during akinesis 25 . 

The mean velocity in the vertical direction was not significantly 
greater than 0 ( — 0.11 ± 0.46, right panel of Fig. 1C), indicating that 
plankton neither emerged from the water surface nor sank to the 
bottom. In the horizontal direction, the mean velocity exhibited a 
slight preference to the right (1.44 ± 0.34, right panel of Fig. 1C) for 
an unknown reason. The mean PSD of the swimming velocity is 
shown in Fig. ID. Log-log plots of the PSD were inversely propor- 
tional to frequency in both the horizontal and vertical directions. The 
slope values of the PSDs calculated by a power-law approximation 
were —0.87 and —0.76 in the horizontal and vertical directions, 
respectively ( — 0.81 on average). 

Behavioural test of medaka. Based on the coordinate data derived 
from swimming D. magna, we generated a raw data model of virtual 
plankton in which the swimming patterns of real plankton were 
faithfully reflected. This model was then used to study the attrac- 
tion and predation behaviour of fish (Fig. 2). The attraction be- 
haviour of medaka is summarised (n = 13) in Fig. 3A, and the 
predation behaviour is summarised in Fig. 3B. In the figures, the 
white bars represent an internal control, which was the mean score 
of the first minute without visual stimuli in each trial. The effects of 
visual stimuli presented during the second minute are shown in 
percentages (the score for the first minute = 100%). In the blank 
control, no visual stimulus was presented during the second minute 



either; no significant increase was observed in the second minute 
compared to the value for the internal control (paired f-test). 

The effects of the raw data model were evaluated by comparisons 
to the blank control. With regard to predation behaviour, the raw 
data model score was significantly and drastically higher than that of 
the blank control (Fig. 3B, unpaired f-test, p = 0.0002, t 26 = -4.284), 
but this was not the case with regard to the attraction behaviour 
(Fig. 3A). These results indicated that the lingering time for fish at 
close range to the cathode ray tube (CRT) display did not signifi- 
cantly increase in response to zooplankton movement, although a 
slight increase was observed. This may be dependent on the visual 
acuity of fish 1719 . Predation behaviour better assesses the effect of 
feeding stimulants as well as the selection method 14,15 . 

We next generated 8 additional types of virtual plankton models in 
keeping with the characteristic parameters found in the swimming 
patterns of D. magna: D. magna, slope, speed, 1// noise, Gaussian 
noise, blue noise, constant speed, and stationary dots models (refer to 
Methods). The frequency distributions for the swimming velocity 
(Fig. IB), mean speed (Fig. 1C), and/or PSD slope value (Fig. ID) 
were used to generate virtual plankton with a circular shape. With 
regard to attraction behaviour (Fig. 3A), none of the scores from the 
virtual plankton models were significantly higher than the blank 
control, again suggesting a lesser potency of attraction behaviour 
when studying feeding behaviour. 

Fig. 3B shows the relative frequency of predation behaviour in 
each experimental group. A comparison of the frequency of preda- 
tion behaviour to that of the blank control using an unpaired f-test 
did not detect significant differences between the Gaussian noise 
(p = 0.117,f 27 = -1.617), blue noise (p = 0.092, t 28 = -1.745), con- 
stant speed (p = 0.055, f 27 = -2.005), and stationary dot (p = 0.243, 
f 2S = —1.193) groups. On the contrary, the D. magna (p < 0.0001, 
t 26 = -5.186), slope (p = 0.0014, t 27 = -3.573), speed (p = 0.0004, 
f 27 = -4.04), and i//-noise (p < 0.0001, f 26 = -5.48) model groups 
exhibited remarkably high frequencies of predation behaviour 
towards the virtual plankton. The scores for predation behaviour 
of these 4 effective virtual plankton were comparable to those from 
the raw data model. 

When the value for the second 1 -minute time bin was compared to 
the value for the internal control of the first 1 -minute time bin, a 
significant increase was observed under every condition except for 
the blank control and the stationary dots (paired f-test). This statist- 
ical analysis indicates whether the moving objects, which suddenly 
appeared on the display, were more or less attractive to the medaka. 
Of all of the visual stimuli, the pink noise-containing models had the 
lowest p-values for predation behaviour (e.g. p = 5.83 X10~ 6 , t I3 = 
-7.320 for the D. magnamodel;p = 8.02X 10~ 6 , t 13 = -7.103 forthe 
i/f-noise model; p = 0.041, t 13 = —2.274 for the Gaussian noise 
model; andp = 0.047, t 14 = —2.174 for the stationary dots). 

Discussion 

Signals with PSD slope values close to —1.0 are termed 1// noise or 
pink noise. The name arises from falling between white noise (1//°) 
and red noise (I If 2 ). White noise is a random signal with a flat power 
spectral density. In other words, the signal contains equal power 
within a fixed bandwidth at any centre frequency. Red noise is com- 
monly known as Brownian noise, which is the type of signal noise 
that is produced by Brownian motion, hence its additional name of 
"random walk noise". Pink noise is detected from numerous signals 
that are generated by various natural and artificial processes, and it 
has been extensively studied in various fields such as electronics 26,27 , 
economics 28 , and psychology 29,30 . In biological fields, heartbeats 31 ' 32 , 
electrocorticograms 3335 , and current flow in ion channels 36,37 exhibit 
pink noise. Although the mean speed of D. magna was in line with 
that reported in recent studies 38 41 , no prior studies have reported 
pink noise in D. magna; a random element to Daphnia swimming 
was previously reported, however 15,40 . 
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Figure 1 | Motion analysis of D. magna. (A) A time series plot of the swimming trajectory of D. magna. In total, 1700 consecutive coordinates (57 s in 
total) of the calculated mass centre of single representative D. magna were plotted in the vertical-horizontal axes. (B) A frequency histogram of the 
swimming velocity derived from 12,288 vectors (256 vectors X 48 D. magna). The swimming velocity was calculated separately in the horizontal (right 
direction as positive value) and vertical (up direction as positive value) axes. (C) Left graph: The mean speed in each direction. Left, right, up-, and 
downward speeds are shown separately. Right graph: The mean velocity in the horizontal and vertical axes. (D) A log-log plot of the mean PSD of the 
swimming velocity in the horizontal (top), vertical (middle), and two-dimensional (bottom) axes. The PSDs were extracted from 256 consecutive vectors 
and averaged (n = 48). Best-fit lines using a power-law approximation are shown in black. 
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Figure 2 | Experimental setup and quantification of medaka behaviour. 

(A) Experimental setup. The test aquarium was attached to a CRT display 
on which virtual plankton were shown. The behaviours of the 3 fish in the 
test aquarium as recorded by a webcam placed above the test aquarium. 

(B) Quantification of medaka behaviour. Motion capture and analysis of 
the medaka in the recorded movies were performed using computer 
software. The centres of the head and mass of fish were determined 
automatically. The angle of the body axis against a line orthogonal to the 
CRT display side was calculated from the 2 centres. An area 7.5 mm in 
width from the CRT display was defined as the prey zone. The residence 
time of the centre of the head in the prey zone was calculated as the measure 
of the attraction behaviour. The time when the head centre was in the prey 
zone and the angle of the body axis was within ±30° was used to calculate 
the residence time as a measure of predation behaviour. 



The common term for the 5 effective types of virtual plankton 
including the raw data model was pink noise. In particular, the II f 
noise model was reconstructed as a pure mathematical model gen- 
erated by the phase-locked loop method 42 , suggesting that pink noise 



is one of the key components for triggering predation behaviour in 
medaka. The asymmetric frequency distribution of the swimming 
velocity found in the vertical direction may not be applicable to 
predation behaviour, however (Fig. IB). 

Pink noise literally includes noise or random elements. To invest- 
igate the effect of such noise and random elements, we examined 
Gaussian (Gaussian white noise) and blue noise models. Blue noise 
has a frequency spectrum such that the power spectral density is 
proportional to the frequency. The scores of predation behaviour 
for the Gaussian and blue noise models were surprisingly not signifi- 
cantly higher than for the blank control (Fig. 3B) suggesting that 
noise elements are not critical for predation behaviour. Stochastic 
resonance is a phenomenon that occurs in a threshold measurement 
system when an appropriate measure of information transfer is max- 
imised in the presence of a non-zero level of stochastic input noise 
thereby lowering the response threshold. The resulting system reso- 
nates at a particular noise level 43 . Mechanisms for determining the 
response threshold of predation behaviour may not be satisfied by 
only stochastic resonance. 

Is predation behaviour tuned to proper specific frequencies within 
pink noise? This hypothesis was supported by a recent study with 
virtual plankton 15 in which bluegill sunfish (Lepomis macrochirus) 
preferred a hop-and-sink motion with a fixed hopping rate at a single 
frequency. Our motion analysis of brine shrimp (Artemia salina), 
which are commonly used as prey for medaka in aquaria, also indi- 
cated that the movement speed was concentrated at a single fre- 
quency (10.74 Hz) in both the horizontal and vertical directions 
(unpublished data). However, the veracity of a hypothesis based 
only on response to specific frequencies was doubtful. Strickler 
et al. reported that a significant preference for the hop-and-sink 
motion was not observed when compared with a liner motion at a 
constant speed 15 . The pink noise effect was notably higher than the 
liner motion model with constant speed (constant speed model), as 
shown in Fig. 3B. Although care should be taken in interpreting data 
from one species of fish, fish may not simply prefer a specific fre- 
quency. These lines of evidence suggest a view that a matrix of vari- 
ables with a wide range of low frequencies, represented by pink noise, 
is important for predation behaviour. 

Random motion has been found throughout different species of 
zooplankton in both fresh water and marine settings 44 . Along with 
other species, the velocity fluctuations in the marine zooplankton 
dinoflagellate (Oxyrrhis marina) appear to indicate pink noise 45 . The 
pink-noise model, derived from the prey-predator interactions 
between D. magna and medaka, may be applicable to prey-predator 
interactions in wide range of fresh water and marine species. In 
addition, pink noise phenomena were also found in the trajectories 
of zebrafish 46 and Escherichia coli 47A>i , suggesting that visual motion 
with pink noise is one characteristic of biological signals. Therefore, 
it is intriguing that predators use biological signals that are ubiqui- 
tously found in living organisms for effective hunting. Further- 
more, the pink-noise idea can be extended to the fishing lure of 
anglerfish or even the prey-predator interactions of terrestrial ani- 
mals. Discrete mathematics, as represented by Fourier analysis or 
fractal analysis 49 , have been successfully used to express several 
important natural phenomena. Prey-predator interactions via pink- 
noise motion are an encouraging field for future research. 

Methods 

Animals. Medaka (Oryzias latipes, black variety 50 ) were purchased from a local pet 
shop (Jumbo Encho, Aichi, Japan). They were housed in a 60 L glass aquaria system 
for at least 7 days before being transferred to the experimental system. The stock 
populations (50 fish per aquarium) were kept in aerated and filtered water at 
26 ± 1°C on a 12/12 h light/dark cycle (lighting on from 08:00-20:00) and were fed 
an artificial dry diet (Tetra Killifish Food, Tetra Japan, Tokyo, Japan). The fish were 
fed at 08:30. The holding water was prepared by mixing de-ionised H 2 0 and artificial 
sea salt (18 g/60 L; Tetra Marine Salt Pro, Tetra Japan, Tokyo, Japan). Adult fish 
(body weight 250 ± 50 mg) were used for all experiments. After the study, the fish 
were transferred to a retirement aquaria and are being maintained for another 
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Figure 3 | Attraction and predation behaviour. Nine types of virtual plankton (raw data, D, magna, slope, speed, 1/f-noise, Gaussian noise, blue noise, 
constant speed, and stationary dot models) were presented to medaka, and the attraction and predation behaviours were calculated. In the blank control, 
no virtual plankton were presented. The data shown represents the mean ± SEM. *p < 0.05, **p < 0.01 (2-tailed unpaired /-test with Bonferroni 
correction) . (A) The mean scores of the attraction behaviour. The values are shown as percentages of the mean score of the first minute, in which no visual 
stimulus was presented. (B) The mean scores of the predation behaviour. The values are shown as percentages of the mean score of the first minute, in 
which no visual stimulus was presented. 



research project. D. magna were kindly donated by Prof. Taisen Iguchi 51 . The donated 
plankton were housed immediately after hatching in a 30 L aquarium {20 L of 
housing water) for 4 days prior to the experiment. The stock populations were 
maintained in the same conditions as the medaka but without aeration and filtration. 
A Chlorella suspension (0.3 ml/ 10 L, Chlorella Industry, Tokyo, Japan) was added to 
the housing tank as a food source once a day at 08:30. This study was approved by the 
committee for Animal Experimentation at the National Institutes of Natural Sciences, 
Japan (approval numbers: 10A044 and 11A019). 

Motion analysis of D. Magna. Cuboid glass aquaria [length of inner side: 9.5 cm 
(horizontal width) X 1.0 cm (horizontal depth) X 9.5 cm (vertical depth)] were used 
as test tanks for the motion analysis of D. Magna. The aquaria were filled with housing 



water (with a vertical depth of 8.0 cm). The temperature of the tank water was 
maintained at 26°C ± 1 °C by air conditioning in the experiment room. The bottom of 
each tank was coated with 1% agarose gel (Nacalai Tesque, Kyoto, Japan) to prevent a 
reflection of the illumination. The lateral sides of the tanks, excluding the side where 
the video camera was positioned, were covered with black rubber to prevent excessive 
illumination. The test tank was placed in a dark room. Illumination at the surface of 
the water was adjusted to 3000 lux using 3 white florescent lamps placed just outside 
each vertical tank surface, excluding the side where the video camera was positioned. 
Ten to 20 D. magna were transferred to the tank, and a digital camera was used to 
record the motion of the plankton (EXILIM EX-F1, Casio Computer, Tokyo, Japan). 
The video images (30 fps, 1280 pixels X 720 pixels) were analysed using DIPP-Motion 
2D software (DITECT, Tokyo, Japan). The coordinates of the centre of mass of the 
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plankton were automatically tracked in each video frame. Data points where plankton 
collided against each other or the tank walls were excluded. To accommodate the 
dimensional limits of the virtual plankton reconstituted on a computer display only 
two-dimensional coordinates (horizontal -vertical axes) were analysed. The 
movement distance per frame was calculated from a series of coordinates and 
converted into the swimming velocity (mm/sec), which was calculated separately for 
the horizontal axis {right direction as positive value) and vertical axis (up direction as 
positive value). The swimming speed was calculated as the absolute value of the 
swimming velocity. A power spectral density (PSD) was extracted from a sequence 
of 256 swimming velocities and calculated using Fast Fourier Transformation (FFT). 
The result is shown as the mean PSD of 48 plankton. The slope of the mean PSD 
(log-log plot) was calculated by a power-law approximation. 

Making the virtual plankton. Virtual planktons were presented on a 15-in CRT 
display with a refresh rate of 70 Hz and a resolution of 1024 X 768 pixels. All stimuli 
were controlled by Psychlops software (C+ + library for developing psychophysical 
visual stimulus, http://psychlops.sourceforge.jp/en/; please refer to our previous 
paper 52 ), running on a Windows PC. Visual stimuli were presented within an area of 
500 X 230 pixels (130 X 60 mm 2 ; interpixel distance: 0.26 mm) located in the centre 
of the display. The coordinate at the bottom-left corner of the area was defined as (0, 
0). The virtual plankton were drawn using 6 white round dots (80 cd'in -2 , Psychlops 
oval function of 3 pixels in diameter) on a black background. The starting coordinates 
of the centres of the dots were (65, 100), (145, 100), (225, 100), (275, 100), (355, 100), 
and (435, 100). The positions of the moving dots were updated every 2 video frames 
according to numerical sequences in data files in which timeline data for the moving 
distance per video frame were described (see below). If the virtual plankton exited the 
area, they re-entered the presentation area from a position diagonal to their exit. The 
positions of the moving dots were reset to the starting position every 10 s. 

To reconstruct virtual plankton on a CRT display, 9 types of numerical sequences 
were prepared using the following procedure. A value less than the interpixel interval 
was added to the next value in the progression. New numerical sequences were 
generated for each trial. The validity of numerical sequences was also examined by 
motion analysis of virtual plankton as described previously for the motion analysis of 
D. magna. 

• 1) Raw data model: The raw velocity data of 71 D. magna were randomly con- 
nected, and 2 progressions (1 each for the horizontal and vertical axes), each of 
which consisted of 26,624 numbers, were prepared. In total, 4096 consecutive 
numbers were randomly clipped from the progressions, and 12 progressions (6 
each for the horizontal and vertical directions) were generated. The binding order 
of the raw data and the 12 progressions were replaced for each trial. 

• 2) D. magna model: Progressions were generated by the Wichmann-Hill algo- 
rithm 53 referencing the discrete probability distributions of the raw velocity data 
from D. magna (refer to Fig. IB). Then, the PSD of each progression was calcu- 
lated by FFT and transformed to fit a function of y — x~ w (N — —0.87 for the 
horizontal axis and —0.76 for the vertical axis) with a convolution filter 54 . The 
mean speed in the horizontal and vertical directions was adjusted to 4.42 and 
5.76 mm/s, respectively. 

• 3) Slope model: Progressions consisting of Gaussian random numbers were 
generated by the Wichmann-Hill algorithm. Then, the PSD of each progression 
was calculated by FFT and transformed to fit a function of y — x~ N (N = —0.87 
for the horizontal axis and —0.76 for the vertical axis) with a convolution filter. 
The mean speed in both the horizontal and vertical directions was adjusted to 
5.09 mm/s. 

• 4) Speed model: Progressions consisting of Gaussian random numbers were 
generated by the Wichmann-Hill algorithm. Next, the PSD of each progression 
was calculated by FFT and transformed to fit a function of y = x~ w (N — — 0.8 1 for 
both axes) with a convolution filter. The mean speed in the horizontal and vertical 
directions was adjusted to 4.42 and 5.76 mm/s, respectively. 

• 5) l//noise: Progressions for which the slopes of the PSDs (power-law approxi- 
mation) equalled —1.0 were generated by the phase-locked loop method 42 . The 
mean speed in both the horizontal and vertical directions was adjusted to 
5.09 mm/s. 

• 6) Gaussian noise: Progressions consisting of Gaussian-distributed random 
numbers were generated by the Wichmann-Hill algorithm. The mean speed in 
both the horizontal and vertical directions was adjusted to 5.09 mm/s. 

• 7) Blue noise: Blue noise sequences were generated using Gaussian-distributed 
random numbers and a convolution filter. The value of N was adjusted to 3.0, and 
the mean speed in both the horizontal and vertical directions was adjusted to 
5.09 mm/s. 

• 8) Constant speed: Each virtual plankton moved with a speed of 5.09 mm/s in a 
constant direction. The direction of motion was decided by a uniform random 
number generator (linear congruential generator 55 ). 

• 9) Stationary dot: All virtual plankton maintained their starting positions during 
the stimulation periods. 

Behavioural test for medaka. The experimental setup is shown in Fig. 2A. Cubic glass 
aquaria (inner side length of 1 5 cm) were used as test tanks in the medaka behavioural 
experiments. To restrict the exterior visual stimulation of the fish, 2 lateral sides were 
covered with white styrene sheets. The test tanks were filled with 1.46 L of medium; 
the water make up was the same as the housing tank. The vertical depth of the water 



was 6.5 cm. The tank water was maintained at a temperature of 26°C ± 1°C by air 
conditioning the experiment room. A total of 142 naive fish were used for ten types of 
visual stimuli (the number of fish in each condition is shown in Fig. 3). 

First, the test tanks were placed in an anterior area of the recording room. Three 
naive fish were randomly selected from stock populations and carefully transferred to 
the test tank (from 09:00-15:00). Twenty-four hours later, the test tank was moved to 
the recording room and attached to the CRT display, and the animals were allowed to 
acclimate for 1 h. The illumination at the bottom of test tank was adjusted to 1500 lux 
with 2 white florescent lamps that were placed just outside the vertical tank surface. 
After the acclimation period, the animals' behaviour was recorded from above using a 
web camera (Qcam-130E; Logicool, Tokyo, Japan) for 2 min. Visual stimuli were 
presented on the CRT display during the last minute. The resulting video images were 
stored in the AVI file format (15 fps, 960 X 720 pixels), and analysed using 
DIPP-Motion 2D software. 

A brief overview of the quantification criteria for the fish behaviour is shown in 
Fig. 2B. The coordinates of the head and mass centres of the fish were automatically 
tracked in each video frame. The head centre was defined as the centre position of the 
2 mass centres of the eyes. The angle of the body axis against a line orthogonal to the 
inner surface of the tank on the CRT display side was calculated from the coordinates 
of head and mass centres. An area 7.5 mm in width from the inner surface of the tank 
on the CRT display side was defined as the prey zone. The width of the area was 
approximately comparable to the length from the front edge to the root of the pectoral 
fins of the fish (mean full length of fish: 25.5 ± 3.4 mm, n = 134). The staying time of 
the centre of the head in the prey zone was calculated as a measure of the attraction 
behaviour. When the head centre was present in the prey zone and the angle of the 
body axis was within ± 30°, the residence time was calculated as a measure of 
predation behaviour. The results of each trial were calculated as the mean values for 3 
fish. The statistical significance of the results was determined using paired (2-tailed) 
or unpaired (2-tailed with Bonferroni correction) f-tests. A cut off of p < 0.05 
was used as the threshold for statistical significance. All values are reported as the 
mean ± SEM. 
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